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Abstract 

We show that the time-dependent particle density n(r,t) and the current density j(r,t) of a 
many-particle system that evolves under the action of external scalar and vector potentials V(f, t) 
and A(f,t) and is initially in the quantum state |-0(O)) , can always be reproduced (under mild 
assumptions) in another many-particle system, with different two-particle interaction, subjected 
to external potentials V'(r, t) and A'(f,t), starting from an initial state |^'(0)), which yields the 
same density and current as \ip(0)). Given the initial state of this other many-particle system, the 
potentials V'(f,t) and A'(f,t) are uniquely determined up to gauge transformations that do not 
alter the initial state. As a special case, we obtain a new and simpler proof of the Runge-Gross 
theorem for time-dependent current density functional theory. This theorem provides a formal basis 
for the application of time-dependent current density functional theory to transport problems. 

PACS numbers: 



I. INTRODUCTION 



Time-dependent density functional theory (TDDFT) attempts to describe the influence 
of many-body interactions on the the time evolution of a quantum many-particle system in 
terms of an effective local potential that depends on a single collective variable, the particle 
density n(r,t). The possibility of such a description is largely based on a fundamental 
theorem, proved by Runge and Gross in 1984, which guarantees the invertibility of 
the mapping from time-dependent potentials V(f, t) to time-dependent densities, for time- 
evolutions that start from a common initial state |-0(O)). More precisely, this theorem 
asserts that two different potentials V(f,t) and V'(f,t) which are both analytic functions 
of t in a neighborhood of the initial time t — 0, and which are not trivially related by 
a gauge transformation (i.e., by the addition of a merely time-dependent constant to the 
scalar potential), cannot give rise to the same density, starting from the same initial state: 
therefore n(f,t) determines V(r,t) uniquely up to a gauge transformation Q|. 

The Runge-Gross theorem has been considerably strengthened a few years ago by the 
proof of another theorem, which we refer to as the van Leeuwen's theorem, which states 
that the time-dependent particle density of a many-particle system that evolves under the 
action of an external potential V(f,t) and is initially in the state | -0 (0) ) , can always be 
reproduced (under mild assumptions) in another many-particle system, with different two- 
particle interaction, subjected to an external potential V'(f, t), starting from an initial state 
\ip'(0)) which yields the same density and divergence of the current density as \ip(0)). Given 
the initial state of this other many-particle system, the potential V'(r,t) is uniquely defined 
up to a gauge transformation. The content of the van Leeuwen theorem is particularly 
interesting in two special cases: (i) if the second many-particle system is non-interacting, 
then the theorem guarantees the possibility of reproducing the time evolution of the density 
of an interacting many-particle system in a non-interacting many-particle system, (ii) if the 
second many-particle system has the same two-particle interaction as the first, then the 
theorem becomes equivalent to the Runge Gross theorem, which is thereby proved in a more 
direct manner. 

Although the TDDFT has been remarkably successful over the 20 years of its history, 
there are several situations in which it seems useful to reformulate the theory in terms of 
the particle current density j(r,t), leading to what is known as time-dependent current 
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density functional theory (TDCDFT) j^]. It has been pointed out that the time- dependent 
exchange-correlation potential, when regarded as a functional of the density, does not admit 
a gradient expansion (this is the so-called "ultra non-locality problem"), whereas the 
exchange- correlation vector potential of TDCDFT does. For this reason there has been 
great interest in applying TDCDFT to situations in which the standard TDDFT is known 
to have problems, such as optical spectra of solids j^EI and dielectric properties of polymer 



chains 



0. 



The reformulation of TDDFT in terms of the current density is particularly 



a, 



attractive to those who wish to apply DFT methods to molecular transport problems 
for the TDCDFT gives direct access to the electrical current density. Finally, it should be 
noted that this theory allows one to dispose of boundary conditions at infinity, such as the 
ones mentioned in footnote Q]. 

Efforts toprovide a rigorous basis for TDCDFT date back to the pioneering work of Ghosh 
and Dhara |9| , where the invertibility (up to a gauge transformation) of the map from scalar 
and vector potentials V(r, t), A(f, t) to particle density and current density n(f, t),j(r, t) was 
first proved under hypotheses similar to those of the original Runge-Gross theorem. The 
Ghosh-Dhara proof of this theorem is considerably more complex than the original proof of 
the Runge-Gross theorem and, furthermore, does not address the issue of representability 
of the interacting current density evolution in a non-interacting system, which is central to 
the foundation of the Kohn-Sham equation. In this paper we overcome both limitations by 
proving the analogue of van Leeuwen's theorem in TDCDFT. The precise statement of the 
theorem is as follows: 

Theorem. Consider a many-particle system described by the time-dependent hamilto- 
nian 

# (*) = E + t] ) 2 + t] ] + E ^ ' (1) 

i L - 1 i<j 

where V(r, t) and A(f, t) are given external scalar and vector potentials, which are analytic 
functions of time in a neighborhood of t = 0, and U is a two-particle interaction. Let 
n(f,t) and j(f,t) be the particle density and current density that evolve under H from a 
given initial state \tp(0)). Then the same density and current density can be obtained from 
another many-particle system, with Hamiltonian 

= E [2^ + £ )) 2 + *)1 + E ^ > ( 2 ) 

i L J i<j 
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starting from an initial state (^'(O)) which yields the same density and current density as 
1-0(0)) at time t — 0. The potentials V'(r,t) and A'(r,t) are uniquely determined by V(f,t), 
A(r,t), \if)(0)), and |^'(0)), up to gauge transformations of the form 

A'(r,t) -> A\r,t)+VA(f,t) , (3) 

where A is an arbitrary regular function of r and t, which satisfies the initial condition 
A(r, 0) = 0. 

The proof of this theorem is presented in the next section, and its physical implications 
are discussed in Section III. 



II. PROOF 

As a first step toward the proof of our theorem let us observe that, given a set of potentials 
V (f, t) , A(r, t) , one can always make a gauge transformation of the form (JHJ) that kills the 
scalar potential at all times. To this end, one simply chooses A(r, t) to be the solution of 
the differential equation 

^-v«t) (4) 

with initial condition A(r, 0) = 0. We will henceforth assume that such a transformation 
has been done in both the unprimed and primed systems so that the scalar potentials V 
and V are zero at all times. 

The current density operator in the unprimed system is given by 



where {A, B} = (AB + BA) denotes the anticommutator of two operators A and B, and 
the velocity operator of the i-th particle is given by 



^ 1 

Vi = — 
m 



Pi + A(ri,t) (6) 

in 

JL 

dfi 

itly on time via the vector potential. A completely analogous expression can be written, 



where pi = —ihSr is the canonical momentum operator. Notice that this depends explic- 
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of course, for the current density of the primed system, the only difference being the re- 
placement of the velocity operator Vi by v[, which is given by Eq. © with A replaced by 
A'. 

Let us denote by j(r,t) and n(f,t) the expectation values of the current density and 
density operators in the quantum state that evolves from the initial state | -0 (0) ) under the 
hamiltonian H. It is not difficult to verify that j(f,t) obeys the equation of motion 



«M = «M + ' <[ H( t ),j( f , t ) ]> 



dt 



dA(f t) r - - i 

n(r, t) K ' ' - j{f, t)x[Vx A(r, £)J + V ■ (a(f, t)) 



(7) 



where 9A Q t '^ and V x A(f, t) are, respectively, the electric and the magnetic field, a is a 
stress tensor operator, defined below, and (...) denotes the average in the quantum state of 
the unprimed system at time t. The stress tensor is explicitly defined as 



cr a p(r,t) = --^{^{^(f-fi)}} - ^^[/(f i ,f i )(5( 



r- ri)8 af} 



(8) 



i j^i 



and its "divergence" is a vector with components [V • a(f, t)] a = ^z£L 



->P drg 



By hypothesis, the current density also obeys the equation 



9j(r ' ^ = n(f, t ) dA Y'^ - j(r, tjxfvx A'(f, t)]+V- (a'(f, t))' , 



(9) 



dt v ' ' dt 

where (...)' denotes the average in the quantum state of the primed system at time t. Notice 
that the stress tensor a' differs from a in two ways: first, because it contains the velocity 
operator v' instead of v and, second, because it contains the two-particle interaction U' 
instead of U. Taking the difference of the two equations (0) and Q we get 

dAA(f, t) 



n(f, t)- 



dt 



j(f,t)x V x AA(r,t) -V • (a'(f,t))' + V • (a(f,t)) 



(10) 



where AA(f,t) = A'(f,t) — A(f,t). Since, by hypothesis, both A(f,t) and A'(r,t) are 
expandable in a Taylor series of time in a neighborhood of t = 0, it follows that their 
difference, AA(f, t) is Taylor-expandable too. We can therefore write 



AA(r,t)=J2AA k (?)t k 



(11) 



k=0 
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where AA k (r) ee & 



t=o 



. Substituting this expansion into Eq. (fTUjl and equating the 



Z-th term of the Taylor expansion on each side of it we easily arrive at 



k=0 



dAA(f, t) 
dt 



k=0 



X 



V x AA 



+ [V-(a(f,t)}'],-[V-(a'(f,t))], 



(12) 



where nfc(r) and j'fc(r) denote the &-th coefficients in the Taylor expansions of n{f) and j(r) 
about t — 0, and, more generally [/(£)]; denotes the Z-th coefficient in the expansion of a 
function f(t) in powers of t about t = 0. The fact that all the quantities appearing in the 
above equation admit such an expansion is a consequence of the analyticity of the vector 
potential and of the time-dependent Schrodinger equation 

M{t)) 



ih- 



dt 



mm)) ■ 



Since 



dAA(f, t) 
dt 



(k + l)AA k+1 {r) 



(13) 



(14) 



we can rewrite Eq. (fT2*|) in the following form: 



l-l i 
n (r)(l + l)AA l+1 (r) = - ^ n,_*(r)(fc + l)AA k+1 (r) + ^ {ji-ft(r) X [VxA4(r]]} 



k=0 



k=0 



+ [V-(a(f,t))'],-[V.(a'(f,t))]z 



(15) 



where we have isolated on the left hand side the k = I term of the sum which originally 
appeared on the left hand side of Eq. (|12j). 

We now show that Eq. (|T3j) is effectively a recursion relation for the coefficients of the 
Taylor expansion of AA(f,t), i.e., a relation that expresses AA;+i(r) in terms of AA k (f), 
with k < I. To this end, we must show that the right hand side of Eq. (JT5|) depends only on 
coefficient AA k (r), with k < I. This is obviously true for the terms in which AA k appears 
explicitly. There are also implicit AA^s, which are hidden in the coefficients of the expansion 
of the density, the current density, and the expectation value of the stress tensor. However, 
the structure of the time-dependent Schrodinger equation (|13|). which is of first order in 
time, guarantees the Z-th coefficient in the Taylor expansion of the quantum states \ifi(t)) 



and \ip'{t)) is entirely determined by coefficients of order k < I in the Taylor expansion of A 
and A': hence all the quantities on the right hand side of Eq. (|15j) are completely determined 
by the coefficients AAk(r), with k < I. (In this argument A is considered a known quantity, 
and A' = A + AA is the quantity we are trying to determine). 

At this point, in order to make the recursion relation Ijl5|) work we only need to determine 
the initial value of AA, namely AAo(f) = A'(r, 0) — A(r, 0). This is easily done, since from 
the equality of the densities and current densities of the primed and unprimed systems it 
follows that 

n(r,0)AA (r) = (<//(0)|J p (r)|^(0)> _ (</,(0)||(r)|^(0)> , (16) 

where j p (r) = ^ J2i ^i) j is the "paramagnetic" current density operator, which 

has the same form in the primed and unprimed system. Thus, the recursion relation (|15jl . 
together with the initial condition completely determines the Taylor expansion of the 
potential A'(r, i) that yields, in the primed system, the same current density that A(f, t) 
yields in the unprimed one. According to our hypotheses, a knowledge of the coefficients of 
the Taylor expansion of A'(r, t) is equivalent to a knowledge of the function A'{r, t) itself. 
This completes the proof of our theorem. 

III. DISCUSSION 

Two special cases of the theorem proved in the previous section deserve a special discus- 
sion. 

1. The primed system coincides with the unprimed system, i.e. U = U' and ^(O)) = 
|^'(0)). In this case Eq. (116(1 above implies that AAq{t) = 0, and then it follows 
from Eq. (USD that AA k (f) = for all k, i.e. A'(f,t) = A(r,t) at all times. This 
important result is just a statement of the Runge-Gross theorem for TDCDFT: it 
asserts that two vector potentials that produce the same current density starting from 
the same initial state of a many-particle system must necessarily coincide, up to a gauge 
transformation: the map from vector potentials to current densities is invertible. As 

n 

noted in the introduction, this theorem was first proved by Ghosh and Dhara [9] by a 
different method, similar in spirit to the original proof of the Runge-Gross theorem. 
The present proof provides a simpler route to the same conclusion. 
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2. The primed system is non interacting, i.e., U' = 0. In this case our theorem provides a 

— * 

resolution of what could be called the non-interacting A-representability problem. In 
other words, the theorem shows that the current density produced by a vector potential 
A in an interacting many-particle system can also be obtained in a non-interacting 
system, under the action of a suitable vector potential A'. This is certainly possible 
if \ip'(0)) = \ip(0)) (in which case we must have A'{r, 0) = A(r, 0)), but it may be 
more generally possible if one chooses for |^'(0)) a single Slater determinant (or even a 
noninteracting ground-state) that yields the correct initial density and current density. 
Thus, the theorem provides a solid basis for the use of the time-dependent Kohn-Sham 
equation, which indeed attempts to reproduce the correct current density in a system 
of noninteracting particles. As pointed out in the introduction, this important result 
lays the foundation for the application of TDCDFT to molecular transport problems. 
Notice, however, that the theorem does not say anything about the possibility of 
producing an arbitrary time-dependent current density by means of a suitable vector 
potential. 
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